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The quantum gravity processes that have taken place in the very early 
Universe are probably responsible for the observed large-scale cosmological 
perturbations. The comparison of the theory with the detected microwave 
background anisotropics favors the conclusion that the very early Universe 
was not driven by a scalar field with whichever scalar field potential. At 
the same time, the observations allow us to conclude that there is a good 
probability of a direct detection of the higher frequency relic gravitational 
waves with the help of the advanced laser interferometers. 

1. Introduction. Fundamental Constants and Units 

The observable quantum gravity effects are usually considered as a subject 
of study for a very remote future. In any case, they are normally expected to 
take place on microscopic, rather than macroscopic, scales. It is remarkable 
that we are probably facing these effects right now, and on extremely large 
scales, in the form of very long- wavelength cosmological perturbations re- 
sponsible for the observed [1] large-angular-scale anisotropics in the cosmic 
microwave background (CMB) radiation. It is also likely that we will have 
direct access to these effects on smaller scales, and in the near future, by 
detecting the relic background of (squeezed) gravitational waves [2] with 
the help of gravity-wave detectors currently under construction or in the 
phase of the design study (for a recent review, see [3]). 
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Since the mentioning of the subject of quantum gravity often instigates 
the feehng of disbehef and panic among the audience, we should first clarify 
our intentions. We are not going to quantize, create and destroy, universes. 
We will be discussing a conceptually simpler problem of quantizing cosmo- 
logical perturbations in the Universe. This approach is quite analogous to, 
say, quantization of excitations in a sample of a condensed matter, or to 
the quantum field description of the squeezed light generation in quantum 
optics experiments [4,5]. 

The cosmological perturbations (density perturbations, rotational per- 
turbations, gravitational waves) consist of a gravitational field component 
and a matter component. The gravitational field and matter variables are 
governed together by a system of coupled differential equations [6] . In case 
of gravitational waves, only the gravitational field component is present. 
An attempt of quantization of the matter part of cosmological perturba- 
tions without quantizing the associated gravitational field part would be 
as inconsistent as an attempt of quantization of the electric components of 
the electromagnetic field without quantizing the magnetic components. 

The central point of our effort is realization of the fact that even if 
there were no cosmological perturbations initially, they must have been 
generated later, under quite general conditions, as a result of parametric 
(superadiabatic) amplification of the zero-point quantum oscillations of the 
perturbation field [7]. The source of amplification, which "pumps" energy 
into the zero-point (vacuum) fluctuations, is the strong variable gravita- 
tional field of the very early Universe. 

Although we are quantizing only the perturbed part of the gravita- 
tional field, we are still quantizing gravity. We expect our results to depend 
explicitly on all the three fundamental constants: G-representing gravity, c- 
representing relativity, and ^-representing quantum theory. Specifically, the 
G and c participate in the definition of the gravitational energy-momentum 
tensor, and the h enters because of our quantum normalization of the en- 
ergy of each mode to ^tko, i.e., to a "half of the quantum in each mode". 
These constants combine naturally in the Planck length Ipi = ^JGhjc? or 
in the Planck mass mpi = yJficjG. The vacuum fluctuations of the fleld in 
question can be visualized as oscillations with a nonzero amplitude propor- 
tional to ^/h. Since everything in this formulation of the problem depend 
on the presence of the vacuum fluctuations, our results must vanish when fi 
is sent to zero (in other words, no quantum theory — no nonzero effects). 
It is important to know why the Planck constant ti participates in our ex- 
pressions, and it is helpful to write and keep track of all the fundamental 
constants. The unfortunate practice of setting outright G = c = h = \ S& 
sometimes an indication that the author is not quite sure why and which 
fundamental constants should be present, and whether, say, the Planck 
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constant should appear in the nominator or in the denominator of the final 
expression. So, for safety, so to say, all the fundamental constants are be- 
ing set to 1. Similarly, the common practice of "measuring" some physical 
quantity (for instance, a scalar field) in units of the Planck mass is not 
at all an element of quantum physics, like the possibility of expressing the 
distance L between the Earth and the Moon in units of the Planck length, 
L ~ 4 X lO'^^/pi, is not an indication of quantum nature of that distance. 

In our study, the fundamental constants G, c, h are present because we 
are quantizing the (perturbed part of) gravitational field and the associated 
matter fields (unless we are considering specifically gravitational waves in 
which case only gravity is being quantized). However, the G, c, h can also 
participate in the final expressions if a discussed problem does not belong 
to the realm of quantum gravity, but involves gravity plus quantum theory. 
A beautiful example is provided by the expression for the maximal masses 
and radii of cold white dwarfs and neutron stars. These objects owe their 
existence to the interplay between gravity and quantum mechanics. It is 
the pressure of the degenerate Fermi gas which through the Pauli exclusion 
principle enables the cold white dwarfs and neutron stars to resist gravi- 
tational collapse [8]. If there were no quantum mechanics {h = 0), these 
objects would not exist. This statement in words is of course supported by 
formulae. The maximal masses of cold white dwarfs (the Chandrasekhar 
mass) and neutron stars are approximately equal to 



XmB J 

where is the mass of a baryon. The equilibrium radius of a white dwarf 
with the mass Mmax is approximately equal to 



where mg is the mass of an electron. These formulae explicitly demonstrate 
that Mmax and R would vanish if the Planck constant were sent to zero. It 
is important to note that the sign of approximate equality in these formulae 
means the approximate equality. It is not a sign of proportionality, which 
would allow many other dimensional quantities to be present, and it is not 
a sign of an estimate based on the dimensional analysis only, which would 
allow the numerical coefficient to be, say, 50 times larger or 50 times smaller 
than 1. A beautiful physics is usually expressed in simple and beautiful 
formulae: the final result depends essentially on the fundamental constants 
and a couple of parameters, such as and me, characterizing the system. 
As we will see below, our predictions for quantum-mechanically generated 
cosmological perturbations in simple cosmological models do also depend 
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essentially on the fundamental constants, combined in Ipi, and a couple of 
parameters characterizing the model. 

The inspection of positions and powers in which the fundamental con- 
stants enter the claimed rcsTilt may serve as a first test of whether we 
are prepared to believe the result. For instance, when you read that "in- 
flation generates" density perturbations because it "stretches the vacuum 
fluctuations beyond the horizon", phis accompanying words about some- 
thing quantum becoming something classical "upon the horizon crossing", 
this explanation may raise your doubts even before any further reading. 
If everything is so simple, why do we consider vacuum fluctuations of the 
"inflaton", and not the vacuum fluctuations of the electromagnetic field 
which we definitely know to exist? Since the vacuum fluctuations of the 
electromagnetic field are being "stretch beyond the horizon" in exactly the 
same manner, we should apparently be able to generate enormous amount 
of photons as well. However, the result for photons is known: strict zero, 
independently on whether the expansion is inflationary or not. This ex- 
ample simply shows that it is the dynamical property of the system, that 
is, whether and how the vacuum fluctuations are coupled to the pumping 
field, and not the universal kinematical property of "stretching" the wave- 
length "beyond the horizon", which is really important. As a first test of 
the inflationary predictions you may wish to check whether and how the 
Planck constant enters the claimed result. Surprisingly, you will not find the 
Planck constant, or when it is present in the form of the Planck mass, it ap- 
pears in the denominator, and not in the nominator, of the final expression. 
Taken for the face value, this position of the Planck constant would sug- 
gest that the effect, owing its existence to the quantum fiuctuations, goes 
to infinity when the Planck constant is sent to zero. Of course, the correct 
position of the fundamental constants is not a guarantee that the result is 
correct numerically. And the lack, or a strange position, of the fundamental 
constants in the final expression is not an indication that the result is nec- 
essarily wrong. The author could have set G = c = ^ = 1 everywhere, or 
could have done this in one part of the participating equations but forgot to 
do the same in the other. However, an unclear physical explanation of the 
expected effect and a chaotic placing of the fundamental constants in the 
resulting expression does raise a suspicion that the numerical result may 
also be wrong. We will return to this point later. 

2. Quantum Normalization 

In order to quantize cosmological perturbations in a traditional field-theoretical 
manner, it is convenient to use the entire general relativity formulated in 
a traditional field-theoretical manner (see [9] and the early papers quoted 



5 



there). The field-theoretical formulation of the Einstein gravity operates 
with the gravitational field potentials (dimensionless, symmetric, sec- 
ond rank tensor) and matter field variables. Both sets of variables are spec- 
ified in the Minkowski space-time 

da'^ = 77^^ dx" dx"" = c^dt^ - dx^^ - dx'^^ - dx^^ . (1) 

The total Lagrangian L of the theory consists of the gravitational part L^, 
which contains only gravitational variables, and the matter part L"*, which 
contains matter variables and gravitational variables. The concrete form of 
and L'^, which we do not write here, distinguishes the Einstein grav- 
ity from many possible alternative theories. The exact Einstein equations 
follow from the variational principle and have the form 

^^,u^^,a _ _ ^ua^^,^. ^ ^^^^ ^ ^^^^t^u ^ ^t^u-^ (3) 

where t^^ is the energy-momentum tensor of the gravitational field, t^^ 
is the energy-momentum of the matter fields including their interaction 
with gravity, and n = SttG/c^. The energy-momentum tensors f^" and r'*'' 
are defined as variational derivatives of and L"^, correspondingly, with 
respect to the Minkowski metric. 

If you want to arrive at the geometrical formulation of the general rel- 
ativity, you are welcome to introduce new functions g^i^ according to the 
rule 

^ g^'^ = T]''^ + h^"' (3) 

and to interpret these new functions as the metric tensor of the curved 
space-time 

ds'^ = dx^ dx" . (4) 

The substitution (3) translates the Lagrangian into the Hilbert-Einstein 

Lagrangian (up to a total derivative) and the field equations (2) into the 
usual geometrical Einstein equations. How the universal gravitational field 
hfj,i, affects equations of motion of matter fields, and how and why, on the 
grounds of physical measurements, we conclude that we "live" in the curved 
space-time (4), is a subject of a separate discussion [9]. 

In addition to the coordinate freedom, i.e., the freedom of choosing 
arbitrary (curvilinear) coordinates in the Minkowski space-time (1), the 
field-theoretical formulation of gravity allows also the gauge freedom, i.e., 
special re-definitions of gravitational and matter variables which do not 
change the field equations (2). We may use this gauge freedom, and we will 
sometimes do that, in order to impose the conditions 



h^"',^, = 



(5) 
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and to further simplify equations (2). 

Prom the field-theoretical perspective, cosmological solutions are noth- 
ing more than specific distributions of gravitational and matter fields in the 
Minkowski space-time (1). Concretely, the homogeneous isotropic FLRW 
models 

ds'^ = c^dt^ - a\t){dx^^ + dx^^ + dx^^) (6) 
are described by the gravitational potentials 

^00 = a^it) - 1 , /ill = /i22 = /i33 = 1 - a{t) . (7) 

Of course, while operating with the Minkowski metric (1), wc do not mean 
that the measurable distance between two objects located at x and x + dx 
is equal to dx. In fact, the gravitational field (7) makes the measurable 
distance equal to a{t)dx. This distance would have been equal to dx if there 
were no gravitational field (7), i.e., for a{t) = 1. But this is a different issue. 
At the moment, we are interested in other sides of the theory. 

In the geometrical approach, the gravitational part of cosmological per- 
turbations is conveniently written as 

ds^ = c^dt^ - a^{t){6ij + hij)dx' dx^ . (8) 

In the field-theoretical approach, we write 

hfiu = + 1 (9) 

where h^^} are given by Eq. (7). The functions hij{ ) participat- 

ing in Eq. (8) and the functions h^i}}{x^,x^,x'^,x^) participating in Eq. (9) 
are totally different objects. They have differing transformation properties 
and differing meanings. But, as functions of {x^ ^), they are simply 

related in the linear approximation that we consider. 

In case of gravitational waves, we can write a spatial Fourier component 
of the field h^^u for a given wave-vector n in the form: 



s=l 



On „— in-x 



(10) 



When writing Eq. (10), the gauge conditions (5) and the remaining gauge 
freedom preserving (5) have been used in order to eliminate h^J and h^^\ 
and to satisfy the requirements 

Pijn^ = Q , PijS'' = Q , (11) 
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where i^ij (s 



1, 2) are constant polarization tensors obeying the conditions 



Pij P 



25. 



Pi 



n 



-n) 



(12) 



In case of gravitational waves, the functions hij and are related as 
h^l^ = a{t)hij, that is, Eq. (10) can also be written as 



hij = - pn(x°) e-^ + Al(a;0) e 



(13) 



After the introduction of a new time parameter r] according to the re- 
lationship dr] = cdt/a{t), and denoting d/dr] by a prime, the linear approx- 
imation to Eq. (2) produces the familiar equation [7] 
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= 



(14) 



for every mode n and for each polarization state s. This equation has two 
linearly independent solutions with two arbitrary numerical coefficients in 
front of them. In the high-frequency regime, that is, when the interaction of 
the waves with the pump field (7) is negligibly small v? S> |^|, the general 
solution to Eq. (14) has the form 



(15) 



Classically, the constants Ai, A2 are arbitrary numbers determined by 
the initial conditions. If there is no classical waves, \Ai\ = \A2\ = 0, there is 
nothing to interact with the pump field, and there is nothing to be amplified. 
Quantum-mechanically, we expect the initial amplitude of the waves to be 
nonzero and to be determined by the zero-point quantum oscillations. We 
want to assign to these words a concrete meaning. We want to determine 
the numerical coefficient C in the general expression for the field operator 
^S]^(^,x): 



c 



(27r)3/2 



2 



Pij(n)-= c„(r/)e'"'^ + c„(r?)e- 



(16) 

where the creation and annihilation operators satisfy the commutation rules 

s s'^ 

[Cn('?), Cj„(r/)] = (5 (n — m)(5ss/, and obey the Heisenberg equations of 
motion with the corresponding Hamiltonian. 

We need to calculate the contribution of gravitational waves to the 
gravitational energy- momentum tensor t^*" (whose exact form is known [9]) 
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in quadratic approximation. Since, at the moment, we are interested only 
in the normahzing coefficient C, we can consider the high-frequency waves. 
Formally, we simply set a(t) = 1 and derive the quadratic approximation 

for t'^'^ using the fact that the functions /iJ^^J satisfy the gauge conditions 
mentioned above, and that the frequency Wn of the wave with the wave 
number n is Wn = cn. Then, the quadratic expression for t'^'^ reduces to 

The quantum-mechanical operator for is constructed by inserting 
Eq. (16) into Eq. (17). The initial vacuum state for each mode of the field 

is defined by the requirement Cn(0)|0) = 0, where Cn(0) and c„(0) are 
the initial values of the annihilation and creation operators, that is, their 
values taken long before any significant interaction with the pump field has 
occurred. The constant C is determined by the requirement that the energy 
of the field is equal initially to a "half of the quantum per mode" , that is 

tood'-^\0) = -n / d^n Wn $^(0|c„(0)S;,(0) + C(0)c„(0)|0) . 

-oo ^ J —oo 

From this requirement we derive C = vT&tt Ipi. Everything in our problem 
is well defined now. 

In case of density perturbations, the gravitational field part and the 
matter part of the perturbations are two sides of the same physical entity. 
For simple models of matter, such as scalar fields and perfect fluids, the full 
set of perturbed Einstein equations can be reduced [10] to a single second 
order differential equation similar to Eq. (14). All the functions describing 
the perturbations can then be found from solutions to this equation via 
simple algebraic and differentiation-integration relationships, to which the 
rest of equations amounts. Concretely, a spatial Fourier component of the 
metric perturbations hij, defined by Eq. (8), can be written in the form 
similar to Eq. (13): 



hij — / ^Pjj 



hniv) e^"-'' + /i„(r?) e 



(18) 



s=l 

s 

where the polarization tensors Pij obeying the conditions (12) are now 

2 

The component hnir]) of the metric perturbations is not an independent 
function, and the perturbed Einstein equations allow us to express this 
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1 

function through the scalar component hniv) of the metric perturbations. 

1 

We denote the scalar component hnir]) = hn{r]) and omit the mode label 
n. 

It is convenient to introduce a new function 7^('?) through the definition 

h+^h'^^ , (19) 

where 



a' 



and H is the Hubble parameter, H = a/a. The function h{ri) changes under 
the small coordinate transformations respecting the form of Eq. (8), but 
the function h + {a/a'^)h' and, hence, the function Jl/a do not change. In 
this sense, the function Jl/a is the coordinate-invariant core of the metric 
component h{rj). 

If the matter is a (minimally coupled) scalar field with arbitrary scalar 
field potential, the single second-order differential equation mentioned above 
is the equation for the function /n(ry), where = ]2: 



a^ 



(20) 



This equation defines the entire dynamical content of the problem and is the 
only one equation which needs to be solved [10]. The components hii{r]) and 

2 

/in (77) of the perturbed metric, as well as the scalar field perturbation, are 
all expressible through the solutions of this equation. The function fl/a is 
reminiscent of the electrodynamic potential allowing to express the electric 
and magnetic components of the electromagnetic field through the solu- 
tions for the potential. It follows from Eq. (20) that in the long-wavelength 
regime the dominant solution is /U ~ That is, the function Jl/a is 

approximately constant, in full analogy with the approximate constancy of 
the metric amplitude fi/a for gravitational waves, Eqs. (13), (14). 

After having found a single independent degree of freedom, completely 
characterizing our coupled system of metric and scalar field perturbations, 
we can proceed with the quantization, in the manner similar to that of 
gravitational waves. It is important to note that if the pump field (7) is such 
(or, in other words, if the scale factor a{rj) is such) that the function 7(7?) 
becomes a constant, Eq. (20) reduces to exactly the same form as Eq. (14) 
for gravitational waves. In any case, the quantum-mechanical operator for 
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hij(r],x) can be written in the form similar to Eq. (16): 

(21) 

2 

For the density perturbations, however, the operator Cn(r/) is not truly 

independent and is expressible through Cn(ry). In its turn, the operator 

Cn(r/) is expressible with the help of Eq. (19) through the operator da{ri) 
for the single degree of freedom, mentioned above, 

c„(7?) = ^ y l{ri)dr,{r]) dri , (22) 

which is subjected to quantization: 

Kiv),dUv)] = S\n-in) , d„(0)|0)=0 . (23) 

The operator dn{rj) is the operator- valued version of the function 'p.-^ij])- 
The normalization constant C in Eq. (21) has been found [10] to be equal 
C = \/247r Ipi- Thus, the quantum normalizing coefficient C for the metric 
perturbations in the case of scalar field density perturbations is practically 
the same as in the case of gravitational waves. 

In this presentation, we are working with the spatially-fiat models (6) 
but the quantum-mechanical generating mechanism is of course indepen- 
dent of this assumption. 



3. Metric Perturbations at the Matter-Dominated Stage 

The (operator- valued) metric perturbations /ijj, Eq. (8), at the matter- 
dominated stage are essentially all we need to know. In case of density 
perturbations and gravitational waves, a synchronous coordinate system 
(8) can, in addition, be made the comoving coordinate system, what we as- 
sume to be done. The functions hij enter the calculation of the microwave 
background anisotropics, and these functions allow us also to find, due to 
the perturbed Einstein equations, other quantities, such as the perturbation 
5p/ p in the matter mass density. In the Schrodinger picture, the quantum- 
mechanical evolution of the initial vacuum state |On)|0_n) results in the 
appearance of the two-mode squeezed vacuum quantum states. With them 
are associated the specific statistical properties of the perturbation field 
itself and the microwave background anisotropies caused by the field [10]. 
However, we will first evaluate the characteristic amplitude of the metric 
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perturbations. The characteristic amphtude h{n) is a dimensionless quan- 
tity defined as the dispersion (square root of variance), or in other words, 
the root mean squared amphtude of the field, per logarithmic frequency 
interval. We use here the Heisenberg picture, in which the initial vacuum 
state |0) remains fixed, while the dynamical evolution is carried by the 
time-dependent field operators. In general, the h{n) is defined through the 
equality 

{0\hijh'^\0) = / h\n)— . 

Jo n 

For gravitational waves, it is convenient to describe the characteristic am- 
plitude h{u) as a function of the frequency v measured in Hz. For density 
perturbations, it is more appropriate to retain the original dependence h{n) 
on the corresponding wave-number (spatial scale) n, since density pertur- 
bations at the matter-dominated stage do not oscillate as functions of time. 

The present matter-dominated stage of cosmological evolution a(r/) ~ rf 
was preceded by the radiation-dominated stage a(ry) ~ rj. How the Universe 
behaved at still earlier times, that is, well before the era of primordial nu- 
cleosynthesis, is not known. Following S. Weinberg, this initial stage of evo- 
lution is called the very early Universe. Usually, we describe the unknown 
evolution of the very early Universe by simple power-law scale factors. These 
are quite representative models, because, if necessary, a more complicated 
evolution can be approximated by a series of power-law pieces. For pur- 
poses of numerical evaluations, it is convenient to regard rj, x^, x^, as 
dimensionless coordinates and to make the scale factor 0(77) carrying the 
dimensionality of length. 

A large class of expanding models is described by the scale factors 

a{r,) = lo\r)\'+^ , (24) 

where rj time grows from —00 and the constant parameter /3 is /? < — 1 
at the initial stage of evolution. The constant Iq has the dimensionality of 
length. With this scale factor the function 7(77) becomes a constant: 

The Einstein equations require the effective equation of state for the matter 
driving the evolution (24) to be in the form 

The case /3 = — 2 corresponds to p = —e and to the de Sitter scale factor 
a{T]) = /okl""^- 
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Having specialized the cosmological evolution a(r/) and knowing the dy- 
namical equations for cosmological perturbations, as well as their initial 
values dictated by the quantum normalization, we can find the charac- 
teristic amplitudes h{iy) and h{n) at any given time rj. When calculating 
the h{v) and h{n), the important property of solutions to Eqs. (14), (20) 
is being used: in the long- wavelength regime, the function fi/a satisfying 
Eq. (14), as well as the function Jl/a satisfying Eq. (20), stay practically 
constant. For definiteness, we will give estimates for the present time, but 
it is clear how the functions h{u) and h{n) scale at earlier times. We will 
also ignore the oscillations in the functions h{v) and h{n) which always 
take place for high-frequency and, correspondingly, short-wavelength parts 
of these spectra as a result of squeezing and standing-wave pattern of the 
perturbations. 

In case of gravitational waves, the results arc as follows. For v < ffj, 
where vh ~ 10^® Hz is the present-day Hubble frequency, 



Specifically for the (dc Sitter) model (3 = —2, the function (27) is inde- 
pendent of frequency v, i.e., we obtain a "flat" part of the spectrum. The 
constant in this model gives the value of the (constant) Hubble parameter 
at the de Sitter stage, H = c/lo, and determines the (constant) curvature 
of the de Sitter space-time. 

In the interval ^ < ^m-, where the frequency Vm is defined by the 
time (energy density) of transition from the radiation-dominated to the 
matter-dominated stage, 10~^^ Hz, 



In the interval v-m ^ ^ ^ where the frequency u,. is defined by 
the time (energy density) of transition from the initial stage of expansion 
(24) to the radiation-dominated stage and above which the spectrum falls 
sharply since the waves with frequencies higher than Uc are not affected by 
the amplification process, fc ~ 10^ Hz in the currently discussed models. 



As we see, the scale factors (24) which are power-law dependent on r] 
time generate spectra which are power-law dependent on frequency v [7]. 
The final result depends on the fundamental constants, combined in Ipi, 
and a couple of cosmological parameters, such as lo and /?, characterizing 




(27) 




(28) 




(29) 
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the model. The sending of the Planck constant to zero would of course 
eliminate the entire expression for h{i'). 

In case of density perturbations, the characteristic amplitude h{n) for 
wavelengths longer than, and comparable with, the Hubble radius is given 
by formula similar to Eq. (27): 

, (30) 

where the wavenumber nu corresponds to the wavelength Xjj = 2TTa/nH 
equal to the Hubble radius. The prcscnt-day \h is Xh 2 x 10^^ cm. 
Specifically for the f3 = —2 model, the function (30) is independent of 
n, that is, we have a "fiat" (or Harrison-Zeldovich-Peebles) part of the 
spectrum. (The slope of the spectrum for density perturbations, but not 
the amplitude (see below) , is a correct part of conclusions in the influential 
early papers [11].) 

The similarity of results (27) and (30) is a consequence of the following 
facts [10]: (i) the basic dynamical equations (14) and (20) strictly coincide 
for the scale factors (24), (ii) the quantum normalizing coefficient C of met- 
ric perturbations is practically the same quantity for density perturbations 
and for gravitational waves, and (iii) the evolution of the long-wavelength 
metric perturbations associated with density perturbations or with gravita- 
tional waves is practically the same (namely, no evolution at all, the metric 
amplitudes stay almost constant) all the way from the initial stage and 
up to the matter-dominated stage (i.e., from the time when a given wave 
"leaves" the Hubble radius and up to the time when it "enters" the Hubble 
radius again). The numerical value of the characteristic amplitude at scales 
comparable with the Hubble radius is approximately the same number, 
both for gravitational waves and for density perturbations. Namely, 

h{uH) ~ ^ (31) 



for gravitational waves, and 



Htih) « ^ (32) 



for density perturbations. This is true independently on how close was the 
parameter /3 to the de Sitter value /? = — 2 at the time when the wavelengths 
of our interest were "leaving" the Hubble radius at the initial stage of 
expansion. This conclusion is in a severe conflict with the conclusion of 
inflationary literature. 
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The conclusion of the inflationary literature is that, at scales comparable 
with the present-day Hubble radius, the amplitude of the perturbation in 
the matter density (5p/p is equal to 



5p_ 
P 



tt2 

^j- , (33) 

H 9o 



where (j)o is the unperturbed scalar field, and the right hand side of Eq. (33) 
is supposed to be evaluated at the time when the wavelength of our interest 

was "leaving" the Hubble radius at the inflationary stage. In the most recent 
literature, one can find the often quoted formula, which explicitly includes 
the Planck mass: 

S'nik) = t;^-^ ^ (34) 



757r2 m% V'^ 

where V is the scalar field potential and V = dV/d(j). 

Leaving aside the lack or strange position of the fundamental constants 
in these formulas, let us first explore them in the vicinity of the most 
favorite (de Sitter) inflationary model. Since the denominator in formulas 
(33), (34) goes to zero in the limit of the de Sitter expansion at the initial 
stage, 4>o ^ 0, F' — *■ 0, these formulas predict arbitrarily large amplitudes 
for 5p/ p today. The same conclusion follows for the characteristic metric 
amplitude, because h{n) ~ 5p/ p when the density perturbations "enter" 
the Hubble radius at the matter-dominated stage. By analyzing the line of 
the inflationary argumentation and reinstating the fundamental constants, 
one can find that, instead of Eqs. (30) and (32), the inflationary prediction 
amounts to ^ 

^(«)infl ~ — "T^ ( — ) (35) 
a/7 lo KuhJ 

for small 7. That is, formula (35) predicts larger and larger amplitudes 
h{n) for models with smaller and smaller values of 7 at the initial stage 
of expansion (i.e., for 2 -|- /3 — > 0, see Eq. (25)). Since the long-wavelength 
metric perturbations are almost constant in time, the divergent as 1/^/7 
behaviour (35) must be postulated as the initial condition (at the first 
"horizon-crossing" ) . 

The inflationary claim about density perturbations leads to an incorrect 
conclusion regarding the possible relative contributions of the qiiantum- 
mechanically generated gravitational waves and density perturbations to 
the observed large-angular-scale anisotropics of CMB. It is known [12,13,14] 
that a numerical value of h{v) at v ^ vh or h{n) at n uh translates into 
a numerical estimate for the quadrupole anisotropy ST/T. According to 
Eqs. (31), (32) the contributions of gravitational waves and density pertur- 
bations are of the same order of magnitude. More accurate evaluation gives 
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some numerical preference to gravitational waves [10]. In contrast, the in- 
flationary claim about density perturbations transforms into the so-called 
consistency relations which are usually formulated as the statement that 
the ratio T/S of the gravity- wave contribution (T) to the contribution of 
density perturbations {S) goes to zero, T/S ^ 0, when the spectral index 
of density perturbations approaches the most favorite (Harrison-Zeldovich- 
Peebles) value, that is, in the limit 2 -|- /3 ^ 0, 7 ^ 0. Obviously, there is 
nothing wrong with the (small) finite contribution of gravitational waves in 
this limit: ST/T Ri Ipi/lo- This result for gravitational waves has never been 
a matter of dispute. It is the claimed divergent contribution of density per- 
turbations which is responsible for T/S* 0. The unacceptable situation 
with the inflationary prediction for density perturbations does not become 
less disturbing if it is formulated as T/S" — rather than S/T 00. The 
main issue is of course the claimed enormous disparity between the (T) and 
(S) contributions, not only the divergence of the (5) contribution by itself, 
because both contributions must be small anyway, in order to be consistent 
with the observations. 

It was shown in detail in [10] (and, briefly, here) that the inflationary 
formula for density perturbations does not follow from correct dynamical 
equations plus correct quantum normalization. (This is one of the "things 
everyone should know about inflation".) However, it is more or less clear 
without detailed calculations that the inflationary proposition can not be 
true. The quantum-mechanical (parametric) generating mechanism can be 
described, both in quantum optics and in cosmology, as the depletion of a 
pump quantum into a pair of signal quanta. A weak laser can not generate 
an arbitrarily large number of squeezed photons at the expense of regulat- 
ing a parameter unrelated to the strength of the pumping light. In contrast, 
the inflationary formula (33) suggests that a relatively weak gravitational 
field determined by a (fixed) can produce an arbitrarily large amount of 
density perturbations at the expense of regulating the denominator in that 
formula, which does not change the strength of the pumping gravitational 
field. The proposition expressed by formulas (33), (34) is known as the 
"standard inflationary result" . In the inflationary literature, ever since this 
claim was formulated on the grounds of dubious arguments, it is regularly 
being "confirmed" . Big confusion has been added by misinterpretation and 
abuse of the "gauge-invariant" formalism, which resulted in the (mathe- 
matically incorrect and physically absurd) statement aimed at justification 
of the divergent "standard result" , namely, that the long-wavelength scalar 
metric perturbations have experienced (in contrast to gravitational waves) 
an arbitrarily "big amplification", proportional to 1/^/7, during a short 
"reheating" transition from p = —e to p = e/3. 
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4. Some Implications of Microwave Background Anisotropics for 
Theory and Experiment 

The observed microwave background anisotropies signify the existence in 
the Universe of small cosmological perturbations with extremely long wave- 
lengths. The wavelengths are of the order of, and longer than, the present- 
day Hubble radius Ih- It is argued [2] that it is difficult to explain the origin 
of these perturbations without invoking the quantum-mechanical (paramet- 
ric) mechanism for their generation. At any rate, in contrast to other pos- 
sibilities, the quantum-mechanical (in fact, quantum gravitational) gener- 
ation is, in a sense, an unavoidable process, since its theoretical foundation 
relies only on general relativity and basic principles of quantum field the- 
ory. This mechanism seems to be sufficient for explanation of the presently 
available data. Moreover, the comparison of theoretical predictions with 
the observations enables us to make certain further conclusions regarding 
the very early Universe and the expected outcomes of new forthcoming 
observations. 

It is important to recall that the mere detection of the quadrupole 
anisotropy at the level ST/T = 5 x 10"^ [1] allows us to conclude that, 
independently on the origin and nature of the responsible perturbations, 
the Universe remains to be homogeneous and isotropic at scales much larger 
than Ih and up to distances about 500 times longer than Ih [13]. At still 
longer scales, the homogeneity and isotropy of the Universe cannot be guar- 
anteed, in the sense that some dimensionless deviations can be larger than 
1 without conflicting the CMB observations [13]. The transition from spa- 
tially fiat cosmological models to open models does not affect this conclusion 
considerably [15]. Also, the mere existence of the long- wavelength cosmo- 
logical perturbations requires them to have special phases and to exist at 
the previous radiation-dominated stage in the form of standing, rather than 
traveling, waves [2]. This conclusion follows from the Einstein equations, if 
we want to propagate the observed perturbations back in time up to, at 
least, the era of primordial nucleosynthesis without destroying the homo- 
geneity and isotropy of that era. The distribution of phases can be only 
very narrow (highly squeezed) with two peaks separated by vr. This distri- 
bution of phases arises inevitably during the generation of squeezed vacuum 
states, and is known in formal quantum mechanics and quantum optics as 
the "phase bifurcation" [16]. 

The processing of the COBE data has resulted in evaluation of the 
power-law spectral index of the long-wavelength perturbations. In the con- 
text of density perturbations, the spectral index is often denoted by n 
(not to be confused with the wavenumber n), and the Harrison-Zeldovich- 
Peebles spectrum (or "flat" spectrum) corresponds to ra = 1. This con- 
vention is caused by the description of the variance for Sp/p in terms of 
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the wavelength intervals, not in terms of the logarithmic wavelength inter- 
vals, used above. The transition to the characteristic metric amplitudes and 
logarithmic intervals would make the Harrison-Zeldovich-Pccblcs spectral 
index equal to zero, in agreement with the "flat" nature of that spectrum. 
In any case, the exact relationship between this spectral index n and the 
parameter /? used above is 

ra = 2/? + 5 , (36) 

so that n = 1 corresponds to /3 = —2. The authors of Ref. [17] derived from 
the COBE data 

n= 1.2 ±0.3 , (37) 

whereas the authors of [18] derived from the same data, but processed in a 
different manner, 

n = 1.84 ± 0.29 . (38) 

In these evaluations, we see the indication that the true value of the 
spectral index n is larger than 1, n > 1. Assuming that the relevant cos- 
mological perturbations were generated quantum-mechanically, the ra > 1 
translates into f3 > —2, see Eq. (36). Specifically, the value n = 1.4, in- 
termediate between the estimates (37) and (38), requires P = —1.8 and 
p = — 1.2e (see Eq. (26)), whereas the conservative value n = 1.2 requires 
P = —1.9 and p = — l.le. With the effective equations of state of this 
kind, p + e < 0, the initial stage of expansion was accompanied by growth 
of energy density and curvature of the space-time. It seems to the author 
that the search for a "microphysical" model of primordial matter capable 
of producing this kind of effective equations of state becomes an important 
theoretical problem. This primordial matter should also be capable of sup- 
porting the scalar-type metric oscillations properly coupled to the pumping 
gravitational field, if we wish to apply the quantum-mechanical generating 
mechanism for the production of density perturbations. It is likely that the 
basic dynamical equation (20) will still be valid, with the replacement of 7 

by ItI- 

Certainly, the required equation of state cannot be accommodated by a 
scalar field, which is the basis of inflationary models. If the spectral index 
n > 1 is confirmed, this will mean that the very early Universe was not 
driven by a scalar field. This conclusion is true irrespective of the form 
of the scalar field potential V{(p). Indeed, the energy density and pressure 
produced by any scalar filed are given by 

^ = ^^' + Vi<l>) > P=2^'^'-m • 
This means that e-|-p=^^^>0, in contrast to the required e + p < 0. 
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The experimentally evaluated spectral index n and, hence, the evalu- 
ated parameter fi enable us to make a relatively firm prediction with regard 
to expected amplitudes of relic (squeezed) gravitational waves in higher 
frequency intervals. This prediction is based on the assumptions that the 
perturbations responsible for 5T/T were generated quantum-mechanically, 
and that the one and the same cosmological model, determined by the de- 
rived fixed /3, is responsible for the entire spectrum, Eqs. (27)-(29). Taking 
into account all known theoretical and observational arguments, it is dif- 
ficult to avoid these input assumptions, but, strictly speaking, they are, 
so far, assumptions. Since the contribution of gravitational waves to the 
large-angular-scale anisotropics is not smaller (in fact, it is somewhat big- 
ger) than the contribution of density perturbations, and since their spectral 
indices are strictly the same (Eqs. (27), (30)), the measured anisotropics 
can be taken as experimental data about gravitational waves. The evahi- 
ated parameter (3 can now be used in formula (29) for /i(z/) in frequency 
intervals accessible for ground-based and space laser interferometers [2]. 

We will start from n = 1.4 and (3 = —1.8. (For some theoretical consid- 
erations about a possibility for n to be close to n = 1.4, see [19].) In this 
case, the predicted signal is well above the expected sensitivity of the pro- 
posed space interferometer LISA [20] and the advanced ground-based LIGO 
(sec [3] and references there). Indeed, the expected signal is h{u) = 10~^^, 
Qg{v) = 10~^ at the LISA-tested frequency v = 10"^ Hz, and h{v) = 10^^^, 
Qglv) = 10-6 at the LIGO/VIRGO/GEO-tested frequency v = 10^ Hz. If 
the more conservative value n = 1.2, (3 = —1.9 is confirmed, this will still 
provide a measurable signal. Namely, h{v) = 10^^*^'^, ^g{y) = 10"^"^ at 
V = 10-3 Hz and h{v) = lO'^^, Q.g{y) = IQ-^O si v = 10^ Hz. 

Thus, in the framework of the theory of quantum-mechanically gener- 
ated cosmological perturbations, gravitational waves are largely responsible 
for the measured large-angular-scale anisotropics, and a more accurate de- 
termination (within certain limits [21]) of the spectral index will make it 
possible to give more accurate estimates for the higher-frequency relic grav- 
itational waves. 

It is necessary to say that the "Pre-Big-Bang" cosmological scenario [22], 
based on a non-Einstein gravitational theory, does also operate with the 
scale factors having (3 > —2. However, according to calculations of the 
authors of this scenario [22], it appears that the quantum-mechanically 
generated gravitational waves exhibit a peak at high frequencies, u = (10^ — 
lO^'^) Hz. The calculated spectrum is far too weak at lower frequencies, and 
is unable to affect and explain the large-angular-scale anisotropics [22]. 

A convincing test of the truly quantum-mechanical origin of certain 
cosmological perturbations should probably exploit the most distinct prop- 
erties of the generated squeezed vacuum quantum states: small variances 
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of phase and large variances of amplitude, and their statistical proper- 
ties in general. The small variances of phase are reflected in the standing- 
wave pattern of the generated fields. This is related to such phenomena as 
the Sakharov oscillations in the density perturbation spectrum and asso- 
ciated oscillations of the higher index multipole components of the CMB 
anisotropies. This property is also responsible for the nonstationary char- 
acter of the relic high-frequency gravitational wave noise, etc. The large 
variances of amplitude are reflected in statistical distributions of various 
measurable quantities, such as the angular correlation variable 
V = f^(ei)f^(e2), etc. 

It appears that we are enforced to deal with the quantum gravity pro- 
cesses, and a detailed comparison of their theoretical predictions with ob- 
servations will certainly be a fascinating area of research in the coming 
years. 
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